OPUCEM: A Library with Error Checking Mechanism for Computing 

Oblique Parameters 

Ozgiir £obanoglu a , Erkcan Ozcan fo , Saleh Sultansoy c and Gokhan Unel rf 

June 29, 2010 



O 

C 

in 

(N 



x 

i 

43 



> 

oo 
(N 

in 

o 
o 



X 



"European Nuclear Research Center, CERN, Geneva, 
Switzerland. ^Department of Physics and Astronomy, 
University College London, London, UK. 'TOBB ETU, 
Physics Department, Ankara, Turkey and Institute of Physics, 
Academy of Sciences, Baku, Azerbaijan. ^Physics and 
Astronomy Department, University of California at Irvine, 
Irvine, USA. 



Abstract 

After a brief review of the electroweak radiative corrections 
to gauge-boson self-energies, otherwise known as the direct 
and oblique corrections, a tool for calculation of the oblique 
parameters is presented. This tool, named OPUCEM, brings 
together formulas from multiple physics models and provides 
an error-checking machinery to improve reliability of numeri- 
cal results. It also sets a novel example for an "open-formula" 
concept, which is an attempt to improve the reliability and 
reproducibility of computations in scientific publications by 
encouraging the authors to open-source their numerical cal- 
culation programs. Finally, we demonstrate the use of OP- 
UCEM in two detailed case studies related to the fourth Stan- 
dard Model family. The first is a generic fourth family study to 
find relations between the parameters compatible with the EW 
precision data and the second is the particular study of the Fla- 
vor Democracy predictions for both Dirac and Majorana-type 
neutrinos. 



1 Introduction 

The categorization of the electroweak (EW) corrections based 
on their contribution types dates back to a study of photon 
propagated four-fermion processes [ ]. The corrections to 
vertices, box diagrams and bremsstrahlung diagrams were all 
considered as "Direct" whereas the propagator corrections due 
to vacuum polarization effects were all named as "Oblique" 
since these participate to the computations in an indirect man- 
ner [2]. 

The EW precision data collected over the last few decades 
by various particle physics experiments have often been used 
to constrain many new models of particle interactions. They 
are particularly useful in checking the allowed parameter 
space of a given model through its contributions especially to 
the vacuum polarization corrections to the boson propagators. 
The main oblique parameters are usually denoted by letters S, 
T, U and the auxiliaries with letters V, W, Y [ ]. As an ex- 
ample, the S parameter estimates the size of the new fermion 
sector and the T parameter measures the isospin symmetry vi- 
olation, i.e. the split between the masses of the new up and 
down-type fermions. The Standard Model is defined by the 
values S = T = U = for a given top quark and Higgs boson 
mass. 

Together with the detailed review by Peskin and Takeuchi 
[2], a number of papers were published, calculating the contri- 



bution of a given model to the oblique parameters. To name a 
few, the estimation of the number of fermion families and neu- 
tral gauge bosons [ ], the validity consideration of the Higgs- 
less models [ ], and the investigation of the Majorana nature of 
the neutrinos [ ] can be cited. However, a number of such pub- 
lications suffer from unusual notations with typos in formulas, 
and errors arising due to utilization of approximations instead 
of exact calculations (with assumptions such as ra# > > mz ) 
or in some cases from unguarded remarks such as "heavy" for 
the new fermions. 

The goal of this work is two fold: the first is to present a 
library to compute the oblique parameters S, T and U both 
with exact one-loop calculations and with some well-defined 
approximations for a number of models and the second is to 
scan the available parameter space for the fourth family mod- 
els. The comparisons between exact and approximate compu- 
tations, and amongst formulas from different papers provide 
an error checking machinery which improves the enduser reli- 
ability. Implemented in C/C++ languages, we call this library 
OPUCEM, which stands for Oblique Parameters Using C with 
Error-checking Machinery. The OPUCEM package consist- 
ing of the library and a set of example driver and presentation 
functions, which are discussed in this manuscript, are publicly 
available [7]. 

The next section describes the technical details of the li- 
brary implementation. Then the following sections are on 
detailed physics studies demonstrating the use cases of OP- 
UCEM. Section 3 uses the OPUCEM library to investigate the 
plausibiliy of a generic fourth Standard Model (SM) family 
(SM4). Section 4 focuses on the implications of the EW data 
from the viewpoint of the flavor democracy (FD) hypothesis, 
which provides a principle theoretical motivation for the po- 
tential existence of a fourth SM family. These case studies 
deal with defining the parameter regions favored by the data 
and define a set of benchmarking points for the SM4. In both 
sections, Dirac and Majorana cases of the fourth SM family 
neutrino are investigated separately. Finally in Section 5, we 
present our concluding remarks. 



2 Function Library Implementation 

The OPUCEM library mainly consists of a header file 
opucem.h declaring the available function prototypes and 
opucem . c implementing them. The functions are grouped by 
the relevant physics cases such as Majorana neutrinos or Higgs 
bosons. In each case, internal comments are used to document 
the source code, indicating the reference paper for each of the 
formula and the nature of the calculations (e.g. exact 1-loop 
calculations, or approximations valid under certain assump- 
tions such as the new fermions being much heavier than the Z- 
boson). Compilation of the library is straightforward, however 
a makefile is provided as is customary. The makefile also fea- 
tures additional targets to produce example command-line and 
graphical-user-interface applications that make use of the li- 
brary in studying the fermions of a fourth SM family [8, 9, 10]. 
One of the primary goals of the implementation is portabil- 
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ity, since we consider portability an auxiliary measure of the 
reliability of the code. To facilitate this goal, the initial im- 
plementation was done in pure C, with the code having no de- 
pendencies or requirements beyond a standards-compliant C 
compiler. Since certain formulas make use of complex num- 
bers, we made use of the complex type found in the C99 stan- 
dard. However, starting with version 00-00-03, the C99 com- 
plex type was dropped in favor of the std::complex template 
in the C++ standard template library. The main motivation 
for this change was to provide an easy roadmap for an even- 
tual replacement of various long double variables with multi- 
precision variables. 

The compilation and execution is compatible with any 
UNIX-like environment with modern versions of the GNU 
compiler collection (gcc) [11], that support a subset of the 
complex functions given in the C99 standard. For example, 
the last version which uses the pure C implementation (vOO- 
00-02) is known to compile with gcc 4.x. On the other hand, 
the current version (v00-00-04), which makes use of the C++ 
std::complex templates, can only be compiled with versions 
of g++ 4.2 onwards, which implements the part of the C++ 
Technical Report 1 (TR1) that pertains to complex numbers. 
TR1 introduced the complex inverse hyperbolic functions of 
C99 into the std::complex template, and these additions are 
expected to appear in the upcoming C++0x standard as well. 

While we do not intend to support Microsoft Windows op- 
erating system as a regular OPUCEM target, we still have 
ported the C version of the library using the GNU develop- 
ment toolchain provided under the Cygwin environment and 
independently by the native tools from the MinGW project 
[12, 13]. The main missing ingredients are complex numbers 
with long double precision. The latest MinGW distribution 
(MinGW 5.1.6) comes with complex.h implementations for 
gcc 3.4.4 and gcc 4.3.4 but only with double precision. For 
this reason, for the Windows port, we rely on an external im- 
plementation of the C99 complex numbers implemented by 
S. Moshier [ ]. It should be noted that this complex-number 
implementation, along with the long double math library from 
the same author, makes it possible to port the C version of OP- 
UCEM to any architecture that supports even rather old ver- 
sions of the GNU Compiler Collection. (For example under 
Cygwin 1.5.25, we were able to compile even with gcc 3.4.4.) 
Under a separate Windows branch, we provide documentation 
and makefiles to facilitate such efforts. 

Finally it should be mentioned that the current version of 
the code, while using the C++ complex templates, is still struc- 
tured in a strictly C-like manner, making no use of the object- 
oriented features of C++. Because of this, it is clear that the 
enduser experience is likely to quickly deteriorate as more 
and more physics models are added to the library in the fu- 
ture. Therefore for the upcoming revisions of the code, the 
planned updates include a better C++ implementation, in par- 
ticular making use of namespaces to separate formulas from 
different physics models. 

The oblique parameter calculation is implemented for the 
following scenarios in the current version of the OPUCEM 
library: 

• New lepton doublets with Dirac [15] or Majorana [16] 
type neutral leptons; 

• New quark doublets with variable 4x4 CKM mixing to 
the third generation [17]; 

• SM and 2 Higgs Doublet Model (2HDM) Higgs bosons 
[18,15]. 

Additionally, a set of functions for recalculation of the SM 
origin in the S — T plane based on reference values of the top 
quark and the Higgs boson masses is also available. 



During the implementation of the library, a number of minor 
errors on the reference papers have been found and corrected. 
These are: 

• mixing angle in the definition of the left-handed weak 
eigenstates of the two Majorana neutrinos in [ 1; their 
equation 10 should read c\ = Mij{M\ + M2) 
Mi/(Mi+M 2 ). 
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• contribution to W and Z boson self energies induced by 
the two Majorana neutrinos in [ ]; the right hand side of 

their equation 3.6 should start as ... = 7yr[[g 
. . . (See the corrected form in [ 1 9]) 



1 n n 2 _ m\-m\ 



2.1 The error-checking mechanism 

To ensure the correct calculation of the oblique parameters, 
a number of validation techniques have been used. The core 
part of the library has been pair programmed, with subsequent 
updates consistently cross-checked against the earlier parts. 
Furthermore, a number of sanity checks have been performed 
during the implementation. Some of these checks are docu- 
mented as embedded comments in the source code. For ex- 
ample, by setting the masses of the second doublet Higgs par- 
ticles in the 2-Higgs-Doublet-Model to zero and taking into 
account the extra multiplicative factor of 3, we have verified 
the calculation of the contribution from the SM Higgs bo- 
son. Other sanity checks have been promoted to become actual 
run-time tests, accessible through clearly indicated safe-to-use 
functions. For example, the contributions to S and U parame- 
ters from Majorana-type neutrinos requires high precision, for 
which the long double type might not be adequate under cer- 
tain values of the input parameters. To catch such cases auto- 
matically, we check whether the computations for Majorana- 
type neutrinos converge to those of the Dirac-type in the fully- 
degenerate limit of the two Majorana masses. Moreover, these 
safe-to-use functions compare approximate and exact com- 
putations and display warnings if absolute or relative differ- 
ence between the two are above certain predefined thresholds. 
Finally, special cases such as equal mass values for up and 
down type fermions, have been addressed carefully with limit- 
ing cases of the formulas derived using symbolic computation 
tools. 

We have also tried to make sure that the OPUCEM imple- 
mentation of the formulas can correctly reproduce the tables 
and plots in the previously published work. An example of this 
has been available in the examples/prd48_pg225 directory 
since the very first public version of the library. The provided 
shell script generates an executable that calculates S and U as a 
function of the lepton mass parameter §, reproducing the first 
figure of [6], as shown in Figure 1. 

2.2 Additional tools and presentation functions 

It is worth noting that some of the comparisons mentioned in 
the previous section cannot go beyond being qualitative anal- 
yses, due to the fact that quantitative information is missing in 
the original publications. In such cases, we have tried to manu- 
ally extract values from the published figures and documented 
the level of agreement as auxiliary notes spread as comments 
in the source code. In a few cases, we have also discovered 
that the representation of data in the figures was incorrect. 

An example of such a mistake is in the Figure 2 of [20], 
where the U — fit results to the EW precision data are being 
plotted on the S — T plane. In that figure, the 90% confidence 
level (CL) ellipse seems to have been plotted instead of the 
2-sigma (95% CL) error ellipse. In our Figure 2, we show 
both the 90% CL and 95% CL error ellipses to highlight the 
enlarged parameter space. In the same paper, we have also 
identified additional minor errors in the value of the T param- 
eter for m//=300 GeV and in the position of the S — T origin, 
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Figure 1: Contributions to S and U due to fourth-generation 
leptons, N\, N2 and £, with masses m\ = % x 50 GeV, 1112 — 
2; x 150 GeV and m t = % x 100 GeV as a function of E, . The top 
plot is from [6] and the bottom plot is the output of OPUCEM 
library. In both cases exact and approximate calculations are 
compared. 



after their recomputation to take into account new reference 
values of the the top quark and presumed Higgs boson masses. 
Although the relative size of all these errors are rather small 
(-3-4%), and the impact on the conclusions of the paper is 
negligible, we mention them here for the sake of accuracy. 

In order to facilitate correct representation of the EW data 
on the S — T plane, we provide our plotting routines as part 
of the OPUCEM package, under the tools/STellipse di- 
rectory. Error ellipses can easily be drawn from the central 
values and the uncertainties of the measured S, T values and 
their correlation coefficient. These plotting functions make use 
of the ROOT package [21], an object-oriented C++ analysis 
framework from CERN and the de facto standard for statisti- 
cal analysis in the experimental high-energy physics commu- 
nity today. 
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Figure 2: Corrected Figure 2 from Kribs et al. [20] 



2.3 Graphical User Interface 

Besides the full-featured standard command-line, the package 
also provides users with a light-weight graphical user interface 
(GUI) for demonstration purposes. The GUI developed within 
the ROOT [ ] framework, is actually an independent appli- 
cation loading the OPUCEM libraries at run-time. It makes 
relevant function calls based on the user set of inputs and sim- 
ply displays the results accordingly on the S-T plane. 

As it can be seen on the Figure 3, the GUI provides sim- 
ple input boxes for the input parameters such as the fermion 
masses and also access to the internal parameters used in the 
calculations such as the sine squared of the weak mixing an- 
gle. The user can perform the S-T calculation, see the resulting 
data points being plotted and save the plot in multiple formats 
(e.g. eps, png, root, etc.) using the buttons at the bottom of the 
GUI window. For example, majority of the S-T plots shown in 
this manuscript have been generated using the OPUCEM GUI. 

Like the library itself, the OPUCEM GUI has been tested 
on various operating systems (OS) running different ROOT re- 
leases. These include but are not limited to SLC-5, Mac OSX 
10.5 and 10.6 running ROOTv5.26.00(b), Ubuntu-9.10 and 
10.4 running ROOTv5. 18.00b. However it should be noted 
that while the GUI can both be interpreted in ROOT and be 
natively compiled, the library functions themselves are always 
loaded dynamically in binary form at run time. 
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Figure 3: The OPUCEM graphical user interface. 

We have chosen the ROOT framework for visualization and 
GUI development as an example, however we do not nece- 
sarily endorse any specific framework. Since OPUCEM is 
a generic function library, it can be loaded and instantly be 
used from within any rapid application development (RAD) 
environment such as MONO[ ] and GAMBAS[ ] for Ba- 
sic, QtCreator[24] for C/C++, and potentially MS Visual Stu- 
dio for .NET pretty much the same way. This makes possible 
applications of the OPUCEM library available on almost any 
platform. 



3 The fourth SM family in General 

For the rest of this paper, we make use of OPUCEM to study 
the Standard Model with 4 generations of fermions (SM4) 
in order to determine the parameter space compatible with 
the precision EW data and to look for possible correlations 
among the SM4 parameters. A fourth SM family, with fermion 
doublets of masses heavier than 100 GeV, could explain some 



emerging new physics hints from current particle physics ex- 
periments and also cure some of the theoretical shortcomings 
of the SM itself [8, 9, 10]. The compatibility of the additional 
generations with the electroweak data was previously investi- 
gated using a multi-parameter fit approach [25, 17]. 

3.1 5, T Formalism and SM4 

Before using OPUCEM to scan the fourth generation parame- 
ter space, we would like to point to the excellent applicability 
of the S, T formalism in studies of SM4, by providing an ex- 
ample comparing the results from a global EW fit (involving 
16 input parameters and more than 10 degrees of freedom) to 
those obtained from the oblique parameters. In a recent pub- 
lication, it has been demonstated that the mixing between the 
third and the fourth generation quarks, 634, plays an important 
role in the determination of the best fit to the EW precision 
data [17]. A resulting plot from this publication is reproduced 
in Figure 4a, in which the most probable value of the Higgs 
boson mass obtained from the global fit is plotted as a func- 
tion of I sin 634 1 . OPUCEM library is able to closely reproduce 
this result without a global fit by scanning the m#-|sin034| 
parameter space, as shown in Figure 4b. The different col- 
ors correspond to different A^"-probability values. There is 
a very good agreement between the most possible values ob- 
tained from OPUCEM and the trend of the global fit results 
from [J 7]. 
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Figure 4: Higgs mass as a function of the mixing angle be- 
tween 3rd and 4th generations, (a) Results of the global fit 
from [17]. (b) A.% 2 probability obtained from a scan of the 
two parameters as obtained using OPUCEM. The 4th gener- 
ation fermion masses are taken to be m e4 = 145 GeV, m V4 = 



3.2 Scanning SM4 with OPUCEM 

Neglecting possible mixings in the fermionic sector as mo- 
tivated by the preference for small | sin 634 1, one deals with 
5 parameters (masses of the four new fermions, M4,o?4,e4, V4, 
and of the Higgs boson) if the neutrino is of Dirac type like 
the other fermions. However in the realm of the SM, there 
is no reason to forbid the Majorana mass terms for right- 
handed components of neutrinos. In this case, the neutrino 
transforms into two Majorana particles. We will denote the 
lighter of these as V4 and the heavier as A/4. Thus in the Majo- 
rana case, one has an extra mass parameter (m/v 4 ) to consider. 
Given the many degrees of freedom, we initially follow a sim- 
ple scanning method: namely, we choose definite value(s) for 
Dirac(Majorana) neutrino mass(es) and scan the masses of the 
three charged fermions and the Higgs boson. We note that 
similar scans with more details have very recently appeared in 
the literature, but only for Dirac-type neutrinos [26, 27]. 
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Figure 5: The correlation between the fourth family quark 
masses and mass differences for Dirac (left) and Majorana 
(right) cases. 

We start by an analysis of the quark sector. The results 
are presented in Figure 5 for Dirac (left) and Majorana (right) 
cases. The mass values, for which the computed S and T fall 
within the 2c error ellipse of the LEP EW working group's 
U = fit [28], are accepted. Darker regions correspond to 
more entries, whereas the white region has no entries. We 
scan the charged fermion masses from 200 GeV to ITeV and 
Higgs mass from 115 GeV to 900 GeV. In the Dirac neutrino 
case, the neutrino mass is assumed to be, m V4 = 90 GeV, the 
lowest value allowed by the LEP data [ ]. It is seen that 
the preferable values for 1*4-0/4 mass difference lie between - 
30 GeV and +40 GeV, independent of the quark mass. How- 
ever, values between -70 GeV and +80 GeV are also allowed. 
Studying other neutrino mass values, we observe that these re- 
sults show a slight dependence on the fourth neutrino Dirac 
mass. In the Majorana neutrino case, we initially assume 
m v . = 80 GeV, the lowest value allowed by the LEP data [ ] 
and OT/y 4 = lOOOGeV. It is seen that although the preferable 
values for U4-CI4 mass difference lie around 100 GeV, other val- 
ues between -HOGeV and +130 GeV are also allowed. If the 
lighter neutrino mass increases, the favoured U4- d.4 mass dif- 
ference also increases, as an example for m V4 = 280 GeV, the 
preferred U4- 0/4 mass difference lies around 170GeV 

In Figure 6, the correlation between the fourth family lep- 
ton mass and quark mass differences are shown for Dirac (left) 
and Majorana (right) neutrinos. In this figure, where the darker 
(lighter) regions correspond to more (less) entries, only the 
mass values falling into the 2a error ellipse of S and T param- 
eters are shown. The charged fermion masses from 200 GeV 
to 1000 GeV and Higgs mass from 115 GeV to 900 GeV have 
been scanned. Neutrino masses are selected as in Figure 5. 

It is seen that only <?4 masses below 275 (390) GeV are 
allowed for Dirac (Majorana) cases for a U4-CI4 mass differ- 
ence that lie between -80 GeV and +90 GeV, independent of 
the quark mass. 
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Figure 6: The correlation between the fourth family lepton 
mass and quark mass differences and for Dirac (left) and Ma- 
jorana (right) cases. 



3.3 Higgs Mass Dependence 

The impact of the Higgs mass is better understood with a re- 
duced set of parameters. To that end, and based on the hint 
from the mass difference studies, the correlation of the mass 
difference between the quark sector and the lepton sector is 
studied for different values of the Higgs boson mass. How- 
ever, in this study, the fermion masses themselves are not 
used as independent variables; while keeping the sum of the 
quark and lepton masses as equal and constant, the probabil- 
ities are calculated as a function of m 4 — c/4 and e\ — v 4 mass 
difference. The upper part of Figure 7 shows the x 2 proba- 
bility of fourth generation's compatibility with the EW pre- 
cision data for ma = 600 GeV but for different values of the 
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fermion mass sums: from left to right m 
600, 1000 and 1600GeV. Different colors correspond to dif- 
ferent probabilities with black being the highest and blue the 
lowest as given in the scale on the right hand side. The 
most probable mass differences constitute an ellipse in the 
m, H — m ( j 4 ,m e4 — m,, 4 plane, independent of the fermion mass 
sums. The ellipse obtained from the fit to the 1600GeV re- 
sults, can be seen on all 3 plots in yellow solid line. The impact 
of the Higgs boson mass is shown in the lower 3 plots, where 
the fermion mass sum is taken as 1600GeV but the Higgs bo- 
son mass changes as 115, 300 and 600GeV, from left to right. 
One can notice that as the Higgs boson mass increases, so do 
the fitted ellipse's radii. 

The results of the fitted ellipse's semi-minor axis length as 
a function of ln(ra//) are given in Figure 8, the upper plot. 
The lower plot shows that the semi-major axis is always 1.75 
times longer than its semi-minor axis. It is therefore possible 
to summarize the results and parameterize the core of the valid 
region as a function of the Higgs mass (all masses in GeV): 
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For Majorana-type neutrinos, the compatibility of the fourth 
generation with EW data depends also on the heavier Majo- 
rana neutrino's mass. This dependency can be illustrated us- 
ing three example scenarios with different Higgs boson mass 
values, shown with open circles in Figure 9. The open cir- 
cle inside the 1-sigma ellipse corresponds to mn =115 GeV 
and m„ 4 =m^ 4 =300GeV, m v , 4 =245GeV and m e4 =355GeV The 
open circle inside the 2-sigma ellipse corresponds to a higher 
Higgs boson mass, of 450 GeV and to fermion masses as 
m„ 4 =335 GeV, m f / 4 =m V4 =265 GeV and m f4 =335 GeV. The third 
open circle, the farthest from the center of the S-T minimum 
is for the heaviest Higgs mass: 900 GeV with the correspond- 
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Figure 7: The correlation between the fourth family lepton 
mass and quark mass differences for Dirac neutrinos. 
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Figure 8: The minor radius and the ratio of the semi -diameters 
for the most probable region ellipses as a function of the Higgs 
mass. 



ing fermion masses as m„ 4 =m f4 =435GeV, m v , 4 =365GeV and 
«Jd 4 =455GeV. For all the three open circles the mass of the 
heavier neutrino, mm, was initially taken to be the same as its 
lighter partner. Keeping their equivalent Dirac mass constant, 
(\f m v 4 x m N 4 = const.), their ratio (myv 4 /m V4 ) was increased 
from 1 to higher values in steps of 0.5 alongside the solid line 
connecting the closed circles. The highest value of the ratio is 
fixed by the lowest experimentally allowed v 4 mass of 80 GeV. 
The plot shows that the increase of the ratio causes an ini- 
tial decrease in T which later becomes an increase. This be- 
haviour is accompanied by a continuous slight decrease in S. 
The combination of these two has the potential of driving the 
test points into the allowed region, even for the m#=900GeV 
case, the farthest from the S-T minimum. It can be concluded 
that the asymmetry between the two Majorana neutrino masses 
provides an opportunity for accommodating very heavy Higgs 
bosons in the Fourth Family model. While the allowed pa- 
rameter space for the Dirac case has recently been consid- 
ered narrow[ ], it is clear that allowing neutrinos to be of 
Majorana-type significantly extends the range of possibilities. 




for various values of the sine of the mixing angle between the 
third and the fourth generations (sin 634). The Higgs boson 
mass was scanned from 115GeV up to 415 GeV in steps of 
50GeV. The ua, mass is set to 610 GeV, compatible with the 
FD predictions and the findings of the previous section. It can 
be seen that a larger mixing angle increases T, compensat- 
ing for larger Higgs boson mass values pulling towards lower 
T and slightly higher S values. The best cases are obtained 
for either small mixing angle (up to | sin 634! =0.05) and small 
Higgs mass (up to 150GeV). Another possibility resulting in 
the 2d error ellipse is attained with a larger mixing angle ( 
I sin 634 J =0.09) and a large Higgs mass (around 300GeV). 



Figure 9: The impact of the Higgs boson and the mass ratio 
of the two Majorana neutrinos. For three example scenarios 
with different Higgs boson masses, the effect of changing the 
neutrino mass ratio is shown on the S-T plane. In each case, 
the open circles represent degenerate Majorana neutrinos (ie. 
the equivalent of a single Dirac-type neutrino). The 1 and 2 a 
error ellipses are from the 2009 results of the U — fit from 
LEPEWWG[37]. 

4 The fourth SM family from Flavor 
Democracy 

Referring the reader to the original [31, 32, 33] and to the re- 
view [34, 35, 36] papers for details, here we summarize the 
main FD results. The FD hypothesis leads the generic case 
of n SM families to n — 1 massless families and one super 
massive SM family with equal fermion masses. Given the 
measured masses of the third SM family fermions, FD re- 
quires the existence of the fourth SM family as the super mas- 
sive one. Mass and mixing patterns of the first three family 
fermions are provided by small deviations from the full fla- 
vor democracy. It should be noted that FD consequences are 
independent on concrete mechanism of fermion's mass gener- 
ation (e.g. Higgs mechanism), the sole requirement is that this 
mechanism should be the same for all SM fermions. 

If the common Yukawa constant of the fermion-Higgs inter- 
actions is taken to be equal to SUw(2) coupling constant giy, 
the fourth family fermion masses are equal to 2\/2gwT] (where 
r/«245 GeV), therefore, m 4 sw 4\/2m w « 450 GeV. On the 
other hand, the fourth family masses are restricted from above 
(OT4 < 1 TeV) due to partial wave unitarity. Depending on the 
nature of the neutrino, again we consider two distinct cases: 

• The fourth family neutrino is of Dirac nature. Then its 
mass should be equal to W4, leading to rrid. ~ m e4 sa 
m V4 m W4. The Higgs boson's mass, and the size of the 
quark mixing between the fourth and the third genera- 
tions are left as free parameters. 

• The fourth family neutrino is of Majorana nature. In 
this case, the two Majorana particles have masses m V4 f=a 
(m4) 2 /M and m^ 4 m M, where M is the Majorana 
mass scale. The current experimental lower bound 



v 4 



> 80GeV leads to M < 2500(6100) GeV if m 4 



450(700) GeV. 

Another possible deviation is due to the large value of the t- 
quark mass: ua, mass is expected to be somehow different than 
niA,. In the remaining of this section we will use OPUCEM 
to investigate the compatibility of these estimations with the 
2009 EW precision data [37]. 

4.1 SM4 with Dirac neutrinos 

For the scenario with Dirac-type neutrinos, we consider an ex- 
ample case with m ( / 4 = m e4 = m V4 = ma, = 550 GeV. The im- 
pact of increasing Higgs boson mass is shown in Figure 10 



ST ellipses 

_ 0.4 

0.3 
0.2 
0.1 



x sin 34 =0.09 
o sin 34 =0.06 
« sin 34 =0.03 
: *sin,,=0.0 



-0.1 
-0.2 
-0.3 




- LEPEWWGSun 

- LEPEWWGSun 



Figure 10: The Higgs mass dependence for various values of 
the third and fourth family quarks mixing angle, | sin 634!. 

A similar approach was taken to study the impact of the 
neutrino mass, m V4 , on the fourth family model's compatibility 
with the EW precision data. The Figure 1 1 contains a scan of 
the Higgs boson mass, from 115 GeV up to 415 GeV in steps 
of 50 GeV for 5 different values of m v4 , 480, 510, 550, 590 and 
620 GeV. The other relevant parameters are m U4 = 610GeV 
and m c [ 4 = m e . = 550 GeV. The third and fourth family quarks 
mixing angle is set as | sin ©34 1 = 0.03 . One can conclude 
that although m V4 = 550 GeV and m# = 1 15 GeV describe the 
parameters most suitable to EW data and also compatible with 
FD, a change in the V4 mass up to ±70 GeV is also allowed as 
long as the Higgs boson mass remains less than 215 GeV. 
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Figure 11: The Higgs mass dependence for various values of 
the fourth family lighter neutrino, V4 . 

The impact of the change in u 4 mass is shown in Figure 12 
for two Higgs boson mass values. The other relevant parame- 
ters are m c ( 4 = m e4 =m V4 =550 GeV and | sin 634! = 0.03. The 
«4 mass was scanned from 480 up to 630 GeV in steps of 
10 GeV. As the 114 mass increases, T value decreases and hits 
a minimum, after that T keeps increasing, meanwhile S de- 
creases slowly but constantly. For larger values of the Higgs 
mass the same behaviour is observed with a starting point con- 
sistent with previous observations, i.e. an increase in S and a 
decrease in T. One can conclude that although the m#=115 
GeV and m L , 4 = 610GeV is most compatible with the EW fits, 
other 114— H mass pairs are also yielding results very close to 



Id ellipse boundary; e.g. mu = 115GeV, m U4 = 490GeV or 
inn = 150GeV, m„ 4 = 610GeV. It is interesting to note that 
the T minimum is obtained for m u . = 1114 = 550 GeV. 
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Figure 14: The 1*4 vs | sin ©34 1 scan. The contributions from 
Higgs bosons with different mass values are shown in different 
colors. 



Figure 12: The fourth family up type quark, 114 mass depen- 
dence for various values of the Higgs mass. 



A cross check was performed by taking the two most com- 
patible M4 mass values found in the previous paragraph and the 
Higgs mass was increased from 115 up to 335 GeV in steps 
of 20 GeV. The other parameters were kept same as before, 
namely m c \ 4 = m e , = m V4 = 550GeV and | sin ©34 1 = 0.03 . 
The results are shown in Figure 13. One can conclude that 
for m U4 =610(490)GeV Higgs masses up to 215(175)GeV 
give results within the 2(7 ellipse. However the smaller Higgs 
masses are favoured. 
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60GeV, independent from m\ mass. For all three m\ values, 
we find that Ax 2 comparable to or lower than the 3-family 
SM's Ax 2 = 1.7 is achievable. 



Table 1: The fourth family parameters at the best Ax 2 values 
compiled for FD hypothesis, Dirac scenario. 



m 4 (GeV) 


400 


550 


700 


m„ 4 (GeV) 


460 


610 


760 


sin 034 


0.03 


0.02 


0.02 


m H (GeV) 


115 


115 


115 


A/ 2 . 

Aftntn 


1.58 


1.67 


1.72 


s 


0.198 


0.202 


0.204 


T 


0.202 


0.202 


0.208 



' ' 



Figure 13: The Higgs boson mass dependence for various val- 
ues of the fourth family up type quark, u\. 



The third and fourth families mixing angle dependence has 
also been investigated. The Figure 14, shows all the mass 
points within the 2(J error ellipse resulting from a scan of m„ 4 
and | sin 034 1 values. The contributions from Higgs bosons 
with different mass values are shown using different mark- 
ers. Other relevant parameters are kept as before, namely 



md 4 



L e^ 



'v 4 



= 550 GeV 



For | sin 634 1 up to 0.07, there are two allowed U4 mass re- 
gions: one above and one below the niA, mass. Inside these 
regions, smaller (larger) Higgs masses necessitate M4 masses 
closer (farther) to the 1114 mass. With the increase of | sin 634! 
value, the two regions merge at 0.07 and continue up to 0.17, 
the maximum allowed value. One should also note that the 
»'» 4 _ 1114 case has a larger allowed range of parameters, it al- 
lows the highest values of the Higgs mass, its larger portion 
survives if the mass points are required to reside in the 1<7 er- 
ror ellipse. Therefore according to FD, and also following the 
observed hierarchy of the second and third generations, one 
expects m„ 4 > m^ . 

Finally for different values of 1114, the free parameters of the 
Dirac case were scanned. The Table 1 contains the mass and 
mixing parameters for the best Ax 2 values obtained from a 
scan of all relevant parameters. It is seen that the favoured 
values are iur = 115 GeV, | sin ©34 1 =0.02 and m u , — nit, = 



4.2 SM4 with Majorana neutrinos 



If the neutrinos are of Majorana type, an additional parame- 
ter, namely the mass of the heavier neutrino {niN 4 ) has to be 
taken into account. This section deals with this case using an 
approach similar to the previous section. The default mass 
pattern is set as: m^. = m e . = m.4 = 550, m V4 = 135, m^ 4 = 
2500 GeV and the Higgs mass greater than 1 15 GeV. In accor- 
dance with the FD predictions, the V4 — N4 mass ratio is se- 
lected to yield a equivalent Dirac mass close to m.4 and the 
U4 — d.4 mass difference is constrained to be less than my/. 
Figure 15 shows a scan of the Higgs boson mass from 150 
up to 900 GeV in steps of 50 GeV for 3 different values of the 
mixing angle between the third and fourth generations. It can 
be seen that a larger mixing angle increases T, compensating 
for larger Higgs boson mass values pulling towards higher S 
and lower T values. The best cases are obtained for either 
small mixing angle (up to | sin 634! =0.05) and Higgs mass be- 
tween 150 and 350 GeV. If the mixing angle is increased up 
to I sin 634 1 =0.1, then a higher Higgs mass range of 550-600 
GeV becomes feasible. 
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Figure 15: The Higgs mass dependence for various values of Figure 17: The fourth family heavier neutrino, N4 mass depen- 
the third and fourth family quarks mixing angle, | sin 634! . dence for various values of the Higgs mass. 



The Figure 1 6 displays the dependency on the V4 and Higgs 
mass values. The Higgs boson mass was scanned from 150 
up to 900 GeV in steps of 50GeV for 4 different V4 masses: 
129, 132, 135 and 138GeV. Other relevant parameters are 
I sin 34 1 = 0.02, m„ 4 =545GeV, m ( , 4 = m e4 = 550GeV and 
OTA/4 = 2500 GeV. One can notice that the effect of increas- 
ing the V4 mass is the opposite of increasing the mixing an- 
gle, I sin 634 1 : it leads to smaller T values. Under the stud- 
ied conditions, the most compatible values are obtained for 
ot V4 = 135 GeV and 150 < m H < 300 GeV. 



Up to this point the u\ — cU mass difference was kept at 
5 GeV. The impact of increasing this difference for different 
values of V4 mass can be seen in Figure 18. The ua, mass 
was scanned from 500 up to 580GeV in steps of 5 GeV for 
three values of V4 mass: 132, 135, 138GeV The other rele- 
vant parameters are ot# = 150 GeV and OTjv 4 = 2800 GeV (one 
of the best pairs found in the previous paragraph), mj 4 — m e4 = 
550GeV and | sin 634 1 = 0.02. The 114 mass increase leads to 
an initial decrease of T and after a minimum it results in an 
increase in T while the S values keep decreasing very slightly. 
For example the minimum for OTv 4 = 135 GeV occurs around 
m U4 = 545 GeV but any value of u\ mass between 530 and 
560 GeV is rather close to the EW data fit minimum. 




Figure 16: The Higgs mass dependence for various values of 
the fourth family lighter neutrino mass, V4 . 
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Figure 18: The fourth family up type quark, u\ mass depen- 
dence for various values of the fourth family lighter neutrino 

mass, V4. 



The Figure 1 7 displays the dependency on the N4 and Higgs 
mass values. The N4 mass was scanned from 1900 up to 
4000 GeV in steps of 100 GeV for four Higgs boson masses: 
115, 150, 200 and 350GeV The other relevant parameters 
are| sin 034 1 = 0.02, ra„ 4 = 545 GeV, ot</ 4 = m e4 = 550GeV and 
ot V4 = 135 GeV. One can notice that the effect of increasing 
the N4 mass is to decrease T sharply and S very mildly. Once 
more the most compatible results (closest to the EW precision 
fit minimum) are obtained for the lightest possible Higgs mass 
(1 15 GeV) for N 4 mass between 2700 and 3300 GeV but also 
the higher Higgs masses (e.g. 350GeV) stays within the lcr 
error ellipse if the N4 mass is between 2100 and 2300 GeV. 



To cross check the optimization of the heavier neutrinos 
mass, otn 4 , it was also scanned from 2500 up to 4000 GeV 
in steps of 100 GeV for two values of U4 mass obtained in 
the previous paragraph. The Figure 19 shows the results of 
this scan. The other relevant parameters are | sin 634 1 = 0.02, 
ot v „ = 135GeV, mj 4 = m e4 = 550GeV and m H = 150GeV It 
can be seen that as in Figure 17, the N4 mass which yields the 
most compatible result is 2800 GeV. 
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Figure 19: The fourth family up type quark, M4 mass depen- 
dence for various values of the fourth family heavier neutrino, 
N A . 

Finally, for different values of m A , all the free parame- 
ters of the Majorana case were scanned with the constraint 
mDimc = yf m N 4 x '^v 4 = '^4 while changing the m^ to m V4 ra- 
tio. The M4 (Higgs) mass was scanned in steps of 5 (10) GeV, 
the mixing angle in steps of 0.01 and the N4 to V4 mass ra- 
tio in steps of 2. The Table 2 contains the mass and mix- 
ing parameters for the best A# 2 values obtained from a scan 
of all relevant parameters. It is seen that the favoured val- 
ues are m# = 115GeV, | sin ©34 1 =0.02, niN 4 /m V4 = 14 and 
m U4 —niA, — 50GeV, independent of yiia, mass. We note that 
the minimum A^ 2 obtained with the scan is much lower than 
the 3-family SM's A^ 2 = 1.7 for all three nu, values. 

Table 2: The fourth family parameters at the best A^ 2 values 
compiled for FD hypothesis, Majorana scenario. 



m 4 (GeV) 


400 


550 


700 


m„ 4 (GeV) 


450 


600 


750 


m N4 (GeW) 


1497 


2058 


2618 


m V4 (GeV) 


107 


147 


187 


| sin 34 1 


0.02 


0.02 


0.02 


m H (GeV) 


115 


115 


115 


Ar 2 . 


0.143 


0.187 


0.215 


s 


0.110 


0.115 


0.117 


T 


0.136 


0.139 


0.146 



5 Conclusions 

We have presented OPUCEM, an open-source function library 
publicly available to calculate the electroweak oblique param- 
eters S, T and U for a number of physics models. This library 
has been heavily tested and shown to be able reproduce var- 
ious earlier results successfully. It provides safe-to-use func- 
tions with an error-checking machinery that relies on compar- 
isons amongst different formulas and between exact calcula- 
tions and well-known approximations. The current version is 
expected to compile, run and give correct results using any 
Unix-like platform with the C++ Technical Report 1 (TR1) 
compiler support (e.g. g++ 4.2 or later). 

OPUCEM is meant to set an example to a new concept that 
we call "open-formula" papers. Many scientific papers include 
numerical examples of various symbolic calculations. With 
time, such examples loose a lot of their value, unless they are 
not readily reproducible by others. Reproduction of old results 
is usually very difficult, since most representations of the nu- 
merical results are on plots and figures, which provide only 
limited precision when one wants to check an independent re- 
implementation of the symbolic calculations. As a solution 
to this, we suggest authors of scientific papers, which present 
and demonstrate new formulas, commit the actual computer 



code of their numerical calculations on a publicly available 
repository. This will also help resolve errors that occasion- 
ally happen when the results are sent to publication, such as 
misprinting of the formulas. 

In addition to the library itself, we have also presented 
an example of its use in studying the mass parameters of a 
fourth generation of Standard Model fermions. Investigations 
yielded regions of the parameter space favored by the precision 
EW data, both for Dirac and Majorana cases of the fourth SM 
family neutrinos. For the Dirac-type neutrinos, the allowed 
mass parameters describe an elliptical ring in the m U4 — m ( i 4 , 
m V4 — m e4 plane when the sum of the fourth generation quark 
masses and the sum of the fourth generation lepton masses are 
kept constant. The radius and the thickness of this elliptical 
ring are dependent on the mass of the Higgs boson and the 
mixing angle between the third and fourth generation quarks. 
We also note that on this ring, EW data favours m ll4 > m c i 4 and 
m V4 < m e4 . For the Majorana case, the additional parameter 
(the mass of the heavier neutrino) gives more flexibility to the 
choice of parameters and especially by adjusting the ratio of 
the masses of the new neutrinos, it is possible to move some 
additional portions of the parameter space to become compati- 
ble with the measured S, T values. Finally, we test predictions 
of the Flavor Democracy hypothesis and find that in either Ma- 
jorana or Dirac case, for the FD predictions to be true, the sine 
of the mixing angle should be about 0.02, m U4 should be about 
50-60 GeV higher than the generic 4th generation mass scale 
and the smaller values of the Higgs mass are preferred. Under 
these conditions, the Aj 2 measured on the S, T plane for an 
FD-motivated 4-generation SM is lower than the 3-generation 
SM. 
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